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We propose a dynamical mean field approach for calculating the electronic structure of strongly 
correlated materials from first principles. The scheme combines the GW method with dynamical 
mean field theory, which enables one to treat strong interaction effects. It avoids the conceptual 
problems inherent to conventional "LDA+DMFT", such as Hubbard interaction parameters and 
double counting terms. We apply a simplified version of the approach to the electronic structure of 
nickel and find encouraging results. 



For systems with moderate Coulomb correlations the 
GW method (and its refinements) jl], [^, |3| is the tool 
of choice for the determination of excited states proper- 
ties from first principles. It is a Green's function-based 
method, in which the effective screened interaction is 
treated at the RPA level, and used to construct an ap- 
proximation to the electronic self-energy. This approach 
cures many of the artifacts encountered when the Kohn- 
Sham orbitals are interpreted as physical excited states, 
while they are actually auxiliary quantities within Den- 
sity Functional Theory (DFT). 

Although the GW approximation (GWA) has provided 
successful treatments of weakly to moderately correlated 
systems such as sp metals and semiconductors, applica- 
tions to more strongly correlated systems with localized 
orbitals indicate a need to go beyond the GWA. For ex- 
ample, in ferromagnetic nickel, it was found Q that the 
GWA is successful at predicting the quasiparticle band- 
narrowing, but does not improve the (too large) exchange 
splitting found in DFT calculations in the local density 
approximation (LDA). The GWA also fails to reproduce 
the 6eV satellite observed in photoemission |2lJ . 

Recently, a new approach to electronic structure cal- 
culations of strongly correlated materials involving d— or 
/— orbitals, has been developed. This approach, dubbed 
"LDA+DMFT" , combines the dynamical mean-field the- 
ory (DMFT) [§ of correlated electron models with DFT- 
LDA calculations Q . It is also a Green's function tech- 
nique, but - unlike GWA - it does not treat the Coulomb 
interaction from first principles. Instead, an effective 
Hamiltonian involving Hubbard-like interaction parame- 
ters in the restricted subset of correlated orbitals is used 
as a starting point. It is thus necessary to introduce 
a "double-counting" correction term. The strength of 
DMFT however is that the onsite electronic interactions 
are treated to all orders, by using a mapping onto a self- 
consistent quantum impurity problem. DMFT has led to 
remarkable advances on electronic structure calculations 
of materials in which the Mott phenomenon or the for- 



mation of local moments play a key role. This is the case, 
e.g. for the satellite structure in Ni, which has recently 
been shown to be correctly described by LDA+DMFT 

The aim of this letter is to take a new step towards 
a first-principles electronic structure calculation method 
for strongly correlated materials. We propose a scheme 
in which the GW treatment of the screened Coulomb 
interaction and exchange self-energy is combined with a 
DMFT calculation for the onsite components of these two 
quantities, in a self-consistent manner ^| . The frequency- 
dependence of the onsite effective interaction (or polar- 
ization) actually requires an extended DMFT scheme (E- 
DMFT), as introduced in earlier work in a model context 
for both the charge and spin channels || |Io|, [yj [L2| . 
This combined GW + (E)DMFT scheme does not make 
use of Hubbard-like interaction parameters and bypasses 
the need for a double-counting correction when imple- 
mented in a self-consistent dynamical manner. In fact, 
using LDA is in principle no longer necessary within such 
a self-consistent implementation. In this work however, 
we implement a simplified version of this scheme on the 
example of ferromagnetic nickel, which serves as a test 
for the feasibility of realistic calculations using this ap- 
proach. 

We consider the Hamiltonian for electrons in a solid in- 
teracting via the Coulomb potential V(r — r') = e 2 /|r — 
r'|. The general strategy of our approach is to con- 
struct a functional of the one-electron Green's func- 
tion G(r,r';r - t') = -(T T ^(r.rtffrV)) and the 
screened Coulomb interaction W Q. Here, T T 
denotes the time ordering operator in imaginary time 
and ^(r,r) [ip^(r, r)] the annihilation [creation] opera- 
tor of an electron at point r at time r. The screened 
Coulomb interaction is defined using the (connected) 
density-density response function: x(r, r';r — r') = 
(T T [p(T,T)-n(r)}[p(v',r')-n(v')]) 

W(r,r';kj n ) = V(r-r')- (1) 
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Following [|13| and [[14) we introduce the free-energy func- 
tional (which generalizes the Luttinger-Ward construc- 
tion) 

T[G,W] = TrlnG - Tr[(Gfr _1 - G _1 )G] 
-^Trln W + ^[{V' 1 - W'^W] + *[G, W] (2) 

In this expression Gh 1 = iuj n + fi + V 2 /2 — v c — Vh is the 
bare Green's function of the solid including the Hartree 
potential Vh(t) = J dr'V(r — r')n(r'). "J is the contri- 
bution to the functional due to electronic correlations 
beyond Hartree. It corresponds to the sum of skeleton 
diagrams which are irreducible with respect to both the 
one-electron propagator and the interaction. 

A formal construction of this functional can be given 
(following [JhJ) by making a Hubbard-Stratonovich trans- 
formation, using auxiliary bosonic fields <p(r, r) conjugate 
to the density fluctuations [^(r, T)tl>(r, r) — n(r)]. The 
effective interaction precisely corresponds to the boson 
correlator: W(r, r', r - t') = (T T 0(r, r)</>(r', t')). The 
functional T is then constructed by a Legendre trans- 
formation with respect to both G and W . A formal 
expression of the correlation functional ^[G, W] (gen- 
eralizing the Luttinger-Ward 3>[G]) can be given, using 
an integration over a coupling constant parameter a be- 
tween the bosonic and fermionic variables: ^[G, W] = 
i / 1 rfa/drdr(< ? !)(r,r)[^ t (r,r)V'(r,r) - n(r)]). The GW 
approximation retains only the first order contribution to 
this functional in the a-expansion, corresponding to the 
exchange diagram % GW = -\TyGWG jfj). 

The equilibrium state of the system corresponds to a 
stationary point of the functional T, which leads to the 
identification of the exchange and correlation self-energy 
Y, xc and of the polarization operator P: 

g = 0^G- = G^-^, E- = g 
ST , , 5$ 
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In the (self-consistent) GW approximation: 
—GW and Pqw = G G (the signs result from the use of 
the Matsubara formalism). 

In order to proceed further, we need to specify a basis 
set. One-particle quantities like G or S are represented 
as G(r,r',^ n ) = Euw^WCW^'^T = 
X^LL'k < / ) L( r )^'ii'(' < , ioJn)4>L' ( r ')* where <fr are localized 
basis functions (e.g. LMTO's) p5[ , centered at an atomic 
position R (and for simplicity assumed to be orthogonal). 
Two-particle quantities such as P or W are represented 
as W(r,r',iu n ) = £ Q/3RR , B*(r)W™' (iv n )B$ (r')*. 
Here B's are linear combinations of cfxp and form an or- 
thonormal set [Q] . Note that the set (jxj) is in general over- 
complete so that the number of B's is smaller or equal 



to the number of </></>. Matrix elements in products of 
LMTOs are then given by 
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with the overlap matrix 0£ = {4>L 1 <pL 2 \B a ) . We note 
that in general we cannot obtain W a p from Wl 1 l 2 l 3 l 4 , 
while the converse is true. 

The functionals T[G, W] and $ [G, W] can thus be 
viewed as functionals of the matrix elements G RR 2 («w n ) 

and W^j 1 (iu> n )- The main idea behind the present work 
is that the dependence of the ^-functional upon the off- 
site components (R ^ R') of G RR and VF RR can be 
treated within the GW approximation, while the depen- 
dence on the onsite components (R = R') requires a 
more accurate treatment. For strongly correlated sys- 
tems, the onsite effective interaction will enter the strong- 
coupling regime in which an RPA treatment is insuffi- 
cient. We thus approximate the functional \& as: 
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(5) 



In this expression, the first term corresponds to the GW- 
functional (written in the specified basis set) and re- 
stricted to off-site components of G and W (i.e associated 
with distinct spheres R ^ R'), namely: 

Gf% ± (t)W?Xl [L > (t)G*:* 2 (-r) (6) 
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GW 



with H^ RR ' L , v given by (Q). Let us note that * 
can also be written as the difference between the com- 
plete GW- functional, and the contributions from the on- 
site components: f™" loc = V GW - ^q"w [G rr , I / F RR ] . 
All the dependence on these onsite components is gath- 
ered into ^i m p- Following (extended) DMFT, this onsite 
part of the functional is generated [ p5[ from a local quan- 
tum impurity problem (defined on a single atomic site), 
with effective action: 

S = J drdr' [-Y, c I( t ^lI(t~t')c l ,(t') (7) 
+ 2 H :c L( T ) c ^W ■■U LiL2L3U (t-t') : c[ 3 (t')c L4 (t') 

where the sums run over all orbital indices L. In this 
expression, is a creation operator associated with or- 
bital L on a given sphere, and the double dots denote 
normal ordering (taking care of Hartree terms). This 
can be viewed as a representability assumption, namely 
that the local components of G and W can be obtained 
from (Q) with suitably chosen values of the auxiliary 
(Weiss) functions Q and U. This is formally analogous 
to the Kohn-Sham representation of the local density 
in a solid. This construction defines the (frequency- 
dependent) Hubbard interactions Ul 1 l 2 l 3 l 11 {u), for a 
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specific material, in a unique manner (for a given ba- 
sis set). Note that Ul 1 l 2 l 3 l 1 must correspond to an 
interaction matrix U a p in the two-particle basis B a via 
a transformation identical to Taking derivatives of 
(|J) as in (|^) it is seen that the complete self-energy and 
polarization operators read: 

T, xc %iu n ) w = ££V(k>„W (8) 

~~ ^^Gw( k ' iuJ n)LL' + [S^ p (iw„)]LL' 
k 

P(q,iU n )afi = P GW {<\,iVn)a(3 (9) 

q 

The meaning of (||) is transparent: the off-site part of the 
self-energy is taken from the GW approximation, whereas 
the onsite part is calculated to all orders from the dy- 
namical impurity model. This treatment thus goes be- 
yond usual E-DMFT, where the lattice self-energy and 
polarization are just taken to be their impurity coun- 
terparts. The second term in (||) substracts the onsite 
component of the GW self-energy thus avoiding double 
counting. As explained below, at self-consistency this 
term can be rewritten as: 

£E£W(t)ll< = - £ W™» l , l ,(t)G l , Li (t) (10) 

k LiL[ 

(where, again, W % ^* L , L , is related to W^ p by an equa- 
tion of the type ^) so that it precisely substracts the 
contribution of the GW diagram to the impurity self- 
energy. Similar considerations apply to the polarization 
operator. 

We now outline the iterative loop which determines 
Q and U self-consistently (and, eventually, the full self- 
energy and polarization operator): 

• The impurity problem (Q) is solved, for a given 
choice of Qll 1 and U a p: the "impurity" Green's 
function GfJ^ p = — (T t cl(t)c^, (t')) s is calculated, 
together with the impurity self-energy T,f mp = 
5^imp/$G lmp = Q- 1 - G^ p . The two-particle 
correlation function xl 1 l 2 l 3 l 4 = (■ 4,! ( r ) c £2 ( T ) :: 
c \ 3 ( t ') c La( t ') : )s must also be evaluated. 

• The impurity effective interaction is constructed as 
follows: 

(ii) 

Here all quantities arc evaluated at the same 
frequency [p6[. The polarization operator of the 
impurity problem is then obtained as: Pimp = 
-2S^ imp /SW imp = U^ 1 - W~ m x p , where the matrix 
inversions are performed in the two-particle basis 
B a . 



• From Eqs. (||) and (JsJ) the full k-dependent 
Green's function G(k, iw n ) and effective interaction 
W(q, iu n ) can be constructed. The self-consistency 
condition is obtained, as in the usual DMFT con- 
text, by requiring that the onsite components of 
these quantities coincide with Gi mp and Wi rnp . In 
practice, this is done by computing the onsite quan- 
tities 

G loc {wn) = £[G H - 1 (k,io; n )-E a;c (k ) ia; n )]- 1 (12) 

k 

Wl oc {Wn) = £[V q - 1 -P(q,^„)]" 1 (13) 
q 

and using them to update the Weiss dynamical 
mean field Q and the impurity model interaction 
U according to: 

Q- 1 = Gy a l + S 4mp (14) 
U- 1 =W{ 1 C +P lmp (15) 

This cycle is iterated until self-consistency for Q and U is 
obtained (as well as on G, W, T, xc and P). Eventually, 
self-consistency over the local electronic density can also 
be implemented, (in a similar way as in LDA+DMFT 
]l6| , p7{ ) by recalculating p(r) from the Green's function 
at the end of the convergence cycle above, and construct- 
ing an updated Hartree potential. This new density is 
used as an input of a new GW calculation, and conver- 
gence over this external loop must be reached. While 
implementing self-consistency within the GWA is known 
to yield unsatisfactory spectra Q , we expect a more fa- 
vorable situation in the proposed GW+DMFT scheme 
since part of the interaction effects are treated to all or- 
ders. 

The practical implementation of the proposed ap- 
proach in a fully dynamical and self-consistent manner 
is an ambitious task, which we regard as a major chal- 
lenge for future research. Here, we only demonstrate the 
feasibility and potential of the approach within a sim- 
plified implementation, which we apply to the electronic 
structure of Nickel. The main simplifications made are: 
(i) The DMFT local treatment is applied only to the d- 
orbitals, (ii) the GW calculation is done only once, in the 
form Q: T*^ w — Glda- W[Glda], from which the non- 
local part of the self-energy is obtained, (iii) we replace 
the dynamical impurity problem by its static limit, solv- 
ing the impurity model (^) for a frequency-independent 
U = U(lo = 0). Instead of the Hartree Hamiltonian 
we start from a one-electron Hamiltonian in the form: 
Hlda - V?™ local - iTr£t mp (0). The non-local part 
of this Hamiltonian coincides with that of the Hartree 
Hamiltonian while its local part is derived from LDA, 
with a double-counting correction of the form proposed 
in in the DMFT context. With this choice the self- 
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FIG. 1: Partial density of states of d-orbitals of nickel (solid 
[dashed] lines give the majority [minority] spin contribution) 
as obtained from the combination of GW and DMFT (see 
text). For comparison with LDA and LDA+DMFT results 
see pi, for experimental spectra see |2l|. 



FIG. 2: Band structure of Ni (minority and majority spin) as 
obtained from the linearization procedure of the GW+DMFT 
self-energy described in the text (dots) in comparison to the 
LDA band structure (dashed lines) and the experimental data 
of |2(J (triangles down) and j2l]] (triangles up). 



consistency condition ( |l2| ) reads: 



non — loc 



(16) 



We have performed finite temperature GW and 
LDA+DMFT calculations (within the LMTO-ASAjl|] 
with 29 irreducible k-points) for ferromagnetic nickel 
(lattice constant 6.654 a.u.), using 4s4p3d4f states, at the 
Matsubara frequencies iu> n corresponding to T — 630K, 
just below the Curie temperature. The resulting self- 
energies are inserted into Eq. (|l6|), which is then used 
to calculate a new Weiss field according to (p3). The 
Green's function G^ c (r) is recalculated from the impu- 
rity effective action by QMC and analytically continued 
using the Maximum Entropy algorithm. The resulting 
spectral function is plotted in Fig. (|lj) . Comparison with 
the LDA+DMFT results in @ shows that the good de- 
scription of the satellite structure, exchange splitting and 
band narrowing is indeed retained within the (simplified) 
GW+DMFT scheme. 

We have also calculated the quasiparticle band struc- 
ture, from the poles of (16), after linearization of 
S(k, iuj n ) around the Fermi level pTj. Fig. (0) shows 
a comparison of GW+DMFT with the LDA and experi- 
mental band structure. It is seen that GW+DMFT cor- 
rectly yields the bandwidth reduction compared to the 
(too large) LDA value and renormalizes the bands in a 
(k-dependent) manner. 

We now discuss further the simplifications made in our 
implementation. Because of the static approximation 
(iii), we could not implement self-consistency on Wi oc 
(Eq. @). We chose the value of U{u = 0) (~ 3.2eV) 
by calculating the correlation function \ an< i ensuring 



that Eq. ([11]) is fulfilled at oj = 0, given the GW value 
for W 1oc (uj = 0) (~ 2.2eV for Nickel Jl|). This pro- 
cedure emphasizes the low-frequency, screened value, of 
the effective interaction. Obviously, the resulting impu- 
rity self-energy Ej mp is then much smaller than the local 
component of the GW self-energy (or than Vj:° c ), espe- 
cially at high frequencies. It is thus essential to choose 
the second term in (^) to be the onsite component of the 
GW self-energy rather than the r.h.s of Eq. (|l0|). For the 
same reason, we included V^° c in Eq.(|l6|) (or, said differ- 
ently, we implemented a mixed scheme which starts from 
the LDA Hamiltonian for the local part, and thus still 
involves a double-counting correction). We expect that 
these limitations can be overcome in a self-consistent im- 
plementation with a frequency-dependent U(lo) (hence 
fulfiling Eq. (|l(])). In practice, it might be sufficient to 
replace the local part of the GW self-energy by S jmp 
for correlated orbitals only. Alternatively, a downfolding 
procedure could be used. 

In conclusion, we have proposed an ab initio dynam- 
ical mean field approach for calculating the electronic 
structure of strongly correlated materials, which com- 
bines GW and DMFT. The scheme aims at avoiding the 
conceptual problems inherent to "LDA+DMFT" meth- 
ods, such as double counting corrections and the use of 
Hubbard parameters assigned to correlated orbitals. A 
full practical implementation of the GW+DMFT scheme 
is a major goal for future research, which requires fur- 
ther work on impurity models with frequency-dependent 
interaction parameters [^2], ^4j as well as studies of 
various possible self-consistency schemes. 

During completion of this work, we learnt about 
Ref. ||] in which a GW correction to the E-DMFT 
scheme has been successfully implemented, in a dynam- 
ical manner, for a one-dimensional extended Hubbard 
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